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Abstract. 

It is observationally as well as theoretically well established that the winds of hot, 
massive OB-stars are highly structured on a broad range of spatial scales. This paper 
first discusses consequences of the small-scale structures associated with the strong 
instability inherent to the line-driving of these winds. We demonstrate the importance 
of a proper treatment of such wind clumping to obtain reliable estimates of mass-loss 
rates, and also show that instability simulations that are perturbed at the lower boundary 
indeed display significant clumping quite close to the wind base, in general agreement 
with observations. 

But a growing subset of massive stars has also been found to possess strong sur- 
face magnetic fields, which may channel the star's outflow and induce also large-scale 
wind structures and cyclic behavior of spectral diagnostics. The paper concludes by 
showing that multi-dimensional, magneto-hydrodynamical wind simulations, together 
with detailed radiative-transfer modeling, can reproduce remarkably well the periodic 
Balmer line emission observed in slowly rotating magnetic O stars like HD 191612. 



1. Introduction 

Hot, massive OB-stars have strong winds, with typical mass-loss rates M « 10~ 6 M Q /yr 
and terminal speeds Voo « 2000km/s, driven by metal-line scattering of the star's in- 
tense continuum radiation field. Th e first quantitative description of such line-driving 
was given in the seminal paper by (ICastor et al.ll 1975L CAK), who assumed a steady- 
state, spherically symmetric, and homogeneous outflow. And since then, extensions 
and applications of this CAK theory have had considerable success in explaining many 
of the observed gross properties and global trends of hot star winds, such as the relation 
between the wind momentum and the stellar luminosity and the basic metallicity de- 
pendence of mass-loss rates. Nevertheless, it has over the years also become very clear 
that, actually, these winds are neither steady-state, nor homogeneous; rather they seem 
to be highly structured on a broad range of sp atial scales. (Fo r a very comprehensive 
review of hot star winds, we refer the reader to iPuls et alll2008h 



This paper summarizes our ongoing efforts in characterizing the nature and con- 
sequences of such wind structure. After first briefly reviewing the background to the 
strong instability inherent to line driving, we next concentrate on confronting theoret- 
ical predictions of the resulting small-scale structure with observations, and examine 
what consequences such wind clumping have for the interpretation of spectral diag- 
nostics. The paper then continues by discussing how the presence of a strong surface 
magnetic field may channel the outflow along closed field lines, and how the large-scale 
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Figure 1. Snapshot of temperature, velocity, and density in a 1-D LDI simulation, 
plotte d versus radius (lef t) and a Lagrangian mass coordinate (right), as defined in 
e.g. lOwocki et al and illustrating how little wind mass is contained in the 

high-speed rarefactions. 



wind structures predicted by simulations of such magnetically confined winds may be 
diagnosed using the rotational phase variation of optical recombination lines like Ha. 



2. Small-scale structure and wind clumping 

The amount of evidence supporting a clumped stellar wind is nowadays overwhelming. 



Sundqvist et al] (1201 lal) give a comprehensive overview of the current status; this paper 
concentrates on two items, namely i) the theoretical background and ii) the effects of 
clumping on wind diagnostics used to derive mass-loss rates and clumping properties. 

2.1. Theory: The line-deshadowing instability 

Linear stability analyses showed already during the 80s ( Owocki & Rybickilll984l[T985 ) 



that the line force responsible for accelerating hot star winds is highly unstable on spa- 
tial scales near and below the Sobolev length L = v t h/(dv/dr) « (v t h/v)r « O.Olr, where 
Vth « lOkm/s is the ion ther mal speed and v » lOOOkm/s is a typical wind flow speed- 
Direc t numerical modeling ( Owocki et all 1 19881 : lFeldmeierlll995k iDessart & Owocki 



120051) has since confirmed that the non-linear growth of this line-deshadowing insta- 
bility (LDI) leads to high speed rarefactions that steepen into strong reverse shocks, 
whereby most wind material is compressed into spatially very narrow "clumps" (or 
shells in 1-D simulations). Fig. Q] illustrates this characteristic structure, which forms 
the basis for our current understanding and interpretation of wind clumping. Further, 
the presence of strong shocks and hot gas in LDI simulations also provide a generally 
accepted expla nation for the soft X- rays observed from OB-stars by orbiting telescopes 
like chandra (Feldmei er et al.|[T997l see also Cohen, this volume). 

The simulation snapshot displayed in Fig. Q] has been calculated using the numer- 
ical hydrodynamics code VH-1 (developed by J. Blondin et al.), with the ra diation line 
force computed following the so-called smooth source function (SSF. lOwocki & Pulsl 
1996) approach. The SSF method allows one to follow the non -linear evol ution of the 



LDI, while also accounting for the stabilizing line-drag effect (lLucvlll984l) of the dif- 
fuse, scattered radiation field. In the displayed simulation, this line drag exactly cancels 
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Figure 2. Left: Density snapshots for an LDI simulation with self-excited struc- 
ture (upper) and one in which the lower boundary has been perturbed by a sound 
wave ( lower). Right: Rad ial stratification of the clumping factor as derived for f 
Pup by Naj arro et al.l (1201 lh vs. theoretical predictions. 



the LDI at the stellar surface, leading to marginal stabilization of the wind base. How- 
ever, as mass parcels move away from the stellar surface, the relative influence of the 
line-drag decreases and structure starts to develop, roughly at r « 1.5/?*. 

But note that the exact amount of such instability-damping in the lower wind de- 
pends critically on delicate details of the competition between the LDI (arising from 
direct absorption of stellar continuum radiation) and the line-drag (arising from the scat- 
tering source function), and so is quite sensitive to, e.g., external perturbations. Indeed, 
we report here on additional model computations in which we perturb the density po at 
the lower boundary by introducing a generic sound wave of amplitude 6p/po = 0.25 and 
period 4000 sec. And such simulations now display significant structure much closer 
to the wind base, as shown by Fig. |2] The right panel of this figure compares the radial 
stratification of the time-averaged clumping factor f c \ = (p 2 )/(p) 2 in simulations with 
self-excited structure and a perturbed lower boundary, respectively, and shows clearly 
how the latter model is more strongly clumped in the inner wind. This result is in 
general good agreement with observations, as further discussed at the end of Sect. 12.21 

2.2. Diagnostics: Deriving mass-loss rates and clumping properties of OB-star 
winds 

The average mass-loss rate from a time-dependent LDI simulation actually is quite sim- 
ilar to that of its steady-state CAK start model, suggesting the resulting wind structure 
only has secondary effect^] on theoretical mass-loss predictions. On the other hand, 
the structure in density and velocity seriously affects the radiative transfer through the 
wind, and thereby also the interpretation of observed spectra. 

To illustrate this, let us compare mass-loss rates derived from different wind di- 
agnostics using smooth wind models. Taking the Galactic 06 supergiant A Cep as an 
example, Fig. [3] shows how the rate derived from the optical recombination line Ha is 
~ 20 times higher than the rate derived from the UV resonance lines of phosphorus v 



'from, e.g., a changed wind ionization balance (Muiir eTet aU20llb . 
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Figure 3. Observed (dotted red lines) and synthetic (bl ack solid and dashed 
lines) PV and Hg l i ne pro files in A Cep. Observations from iFullerton et al. (2006) 
and [Markova et al. ( 2005). Synthetic profiles calculated u sing the unified (pho- 
tosphere+wind) NLTE model atmos phere code fastwind dPuls et alj 120051) and 
the Monte-Carlo code developed by ISundqvist et alj d2010l 1201 lbl) . All profiles 
were calculated assuming a smooth wind, and illustrate how the mass-loss rate 
(6 x 10~ 6 M Q /yr) required to fit Ha is ~ 20 times higher than that required to fit 
the PV lines. 



(PV), the latter also being almost an order of magnitude lower than predicted by the 
line-driven wind theory! Such inconsistencies have been found by many studies (with 
perhaps the most p rominent example being the study of PV in 40 Galactic O stars by 
Fullerton et al. I I2006I) . and are likely a consequence of neglecting wind clumping when 



deriving the mass-loss rates. 

Most diagnostic studies of clumping have assumed that the wind consists of sta- 
tistically distributed, optically thin clumps embedded in a void background medium, 
and neglected any disturbances on the velocity field. The most important result of such 
optically thin clumping is that the opacities of diagnostics depending on the square 
of the density, such as Ha and thermal infra-red/radio emission, are higher than in a 
smooth wind model with the same mass-loss rate. The scaling invariant is yffc\M, 
leading to overestimates of mass-loss rates derived from smooth wind models by -Jfa. 
In contrast, diagnostics depending only linearly on density, such as UV resonance lines 
and bound-free attenuation of X-rays, are not dir ectly affected by optically thin clump - 
ingP[ Combined optical/infra-red/radio studies dPuls et al.l 120061 : [Najarro et all 12011 ) 
using this approach typically show that mass-loss rates derived from Ha and smooth 
wind models should be scaled down by factors 2 ... 5. Such reducti ons have also been 
confirmed by independent studies of bound-free X-ray attenuation dCohen et al.ll2O10L 
201 ll . see also Cohen, this volume), and would imply modest reductions of present-day 
theoretical predictions of mass-loss rates, as well as clumping factors on order f c \ « 10. 



Clumping in the inner wind. While such characteristic clum ping factors are in gen- 
eral accordance with theoretical predictions, diagnostic studies dEversberg et al.|[l998l : 
iBouret et"ai] 120031 : iPuls et all 120061 : ICohen et all l201ll) also suggest the wind is sig- 



2 they can be indirectly affected though, through a modified ionization balance (e.g. lBouret e t al. 2003) 
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Figure 4. Same as Fig. [3] but now accounting for clumping and optically thick 
vorosity (solid lines, see text), and using a single mass-loss rate 1.5 x lO~ 6 M /yr. 
As an illustration of the importance of optically thick clumping for UV resonance 
lines, the dashed lines show synthetic PV profiles calculated with the same mass-loss 
rate but now under the assumption of optically thin clumps (i.e. neglecting vorosity). 



nificantly structured well below the r « 1.5/?* predicted for the onset of clumping by 
standard LDI simulations with self-excited structure. The right panel of Fig.|2]compares 
the radial strat ification of fa derived for £ Pup by the comprehensive multi-diagnostic 
study by Naj arro et all (1201 ll) with theoretical predictions. The figure shows that per- 
turbing the lower boundary (see Sect. 12.11 ) can indeed induce significant wind structure 
also in the lower wind, leading then to quite good overall agreement with observations. 
A journal paper discussing such base perturbations, as well as other effects that may 
affect the predicted structure in the inner wind, is currently in preparation. 



Optically thick clumping - porosity and vorosity. But note that such optically thin 
clumping does not really help explain the unexpected weakness of the PV lines dis- 
played in Fig. HI as these line are basically unaffected by this approach. It is important 
to realize, however, that if the assertion of optically thin clumps is not met for the inves- 
tigated process, additional effects become important in the radiative transfer. For con- 
tinuum diagnostics, self-shielding of opacity within optically thick clumps leads to in- 
creased escape of photons through porous channels in between the clu mps. Such poros- 
ity th u s has the general eff ect of reducing the wind's effective opacity dFeldmeier et al. 



12003k lOwocM et al.l 12004b . Though it has bee n suggested porosity might be impor- 
tant f or X-ray line attenuation in hot star winds dFeldmeier et al.ll2~003l : lOskinova et al 



120061) . recent studies show that thi s is quite unlikely as there is no apparent evidence for 
significant porosity in the data (Cohen et al.ll2008l Leutenegger et al., submitted) and 
it woul d further imply mean free paths between clumps much longer than predicted by 
theory dOwocki & Cohenll2006l : ISundqvist et al.ll2012ah . 

But the situation is quite different for the inherent ly very strong UV re sonance 
lines, in which clumps easily can become optically thick ( Qskinova et al.ll2007T) . Due to 
the narrow Doppler width of the line profile, line photons in a rapidly expanding stellar 
wind can only interact with the plasma over a very narrow spatial range. Therefore 
many such line photons will simply leak through the wind via "porous" channels in 
the velocity field, without ever interacting with the optically thick clumps. Hence this 
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Figure 5. Left: Contour maps of density (upper) and radial velocity (lower) in an 
MHD wind simulation, obtained by aver aging over ~ 100 sna pshots well after initi- 
ation. Right: Observed (ephemeris from lHowarth et al]l2007l) and s ynthetic Ha dy- 
namic spectra for HD 19 162, calculated as described in the text and in Sundqvist et al.l 
J20l2bl) . 



effect has been dubbed velocity porosity, or "vorosity" ( Owockill2008l) . And in analogy 
with spatial porosity, the main effect of such vorosity is to reduce the wind's effective 
opacity. 

Using models accounting for both clumping and optically thick vorosity, Fig. |4] 
again shows PV and Ha fits to A Cep. And in contrast to before (Fig. [3]), both diag- 
nostics can now be well fit using the same mass-loss rate, which now is only a modest 
factor of two lower than the theoretical prediction by lVink et all ( 2000l) . This illustrates 
the importance of a proper treatment of wind clumping when deriving mass-loss rates 
from OB -stars. 



3. Large-scale wind structure due to magnetic fields 

Let us now switch focus from the small-scale structures associated with the LDI to 
the more large-scale wind structures associated with the presence of a strong surface 
magnetic field. Due to recent advances in spectropolarimetric techniques, such large- 
scale fields are being detected in OB -stars at an increasing rate (for an overview of the 
current observational status, see e.g. Wade, this volume). 

3.1. Theory: Dynamical magnetospheres 

The effectiveness of the magnetic field in channeling the stellar wind outflow may be 
characterized by the ratio of the magnetic to wind kinetic energy density: 

B 2 /(Sn) (r/R+y 4 
pv z /2 v/Vco 

where the second equality assumes a dipolar field and defines the so-ca lled wind con- 
finement parameter ?y* = (B^/?^)/(Mvoo) dud-Doula & Owockil |2002| ). with the 
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equatorial surface field strength. A key point is that for 77* > 1, the Alfven radius 
Ra ~ tl+^R-k at which the magnetic and wind kinetic energy densities are equal is lo- 
cated away from the stellar surface, allowing then for wind plasma to be channeled 
along closed field loops toward the magnetic equator. But note also that the steeper 
radial decline of the magnetic energy density than the wind kinetic energy density 
inevitably means the wind will always dominate at large enough radii (as is directly 
evident from the second equality in eqn. [T]) and force the field lines to open up and 
essentially follow the radial flow. 

Magneto-hydrodynamical (MHD) simulations following the time-evolution of the 
competition b etween the magnetic field and the radiatively driven wind confirms this 
basic picture (lud-Doula & Owockll2002l) . The MHD simulation (computed by A. ud- 
Doula) used to synthes ize the Ha emission in the following section is described in 



Sundqvist et al.l (1201 2bi ) and assumes 77* = 50. The left panels of Fig. [5] show that 



below the Alfven radius Ra ~ 2.72?* the wind indeed becomes trapped and the material 
is pulled back onto the star by gravity. But despite this very dynamic behavior, with gas 
constantly in upward and downward motion, the transient suspension of plasma within 
such closed field loops leads to a region around the magnetic equator that statistically 
is overdense. 

Note that this structure is physically dist inct from that predicted for rapidly rotat- 



ing m agnetic stars characterized by Ra > Rk dTownsend & Owockil2005l : lud-Doula et al. 



2008), where R^ is the Kepler corotation radius at which the centrifugal force balances 



gravity. Fur such stars, centrifugal forces can support any trapped material above 7?k, 
which allows for accumulation of wind plasma over very long time-scales and to the 
formation of a centrifugal magnetosphere (CM). In contrast, the characteristic structure 
described above, appropriate for slowly rotating magnetic stars with /?* < Ra < Rk, 
forms a dynamical magnetosphere (DM). An important distinction between a CM and a 
DM is that the latter accumulates wind plasma only over the relatively short dynamical 
time-scale, implying that suc h DMs require the quit e high mass-loss feeding rate of an 
O star in order to emit in Ha (Sundq vist et al.ll2012bl . see also Petit & Owocki, in prep.). 



3.2. Diagnostics: Rotational phase variation of Balmer line emission 

We compute synthetic Ho- profiles directly from the MHD simulation discussed above 
by solving the formal integral of the radiative transfer equation in a 3-D cylindrical 
coordinate system aligned toward the observer. Since we target slowly rotating stars, 
we may neglect rotational effects on the wind dynamics and us e the same simulation for 
any obliquity angle /3 between the rotation and magnetic axes. Sundqvist et al.1 d2012bl) 
give a full description of the computation method. The right panel of Fig. [5] compares 
the rotational phase variation of such synthetic Ha spectra with observ ations of the 
very slowly rotating (P = 538 days, v rot « 1 km/s. iHowarth et al.ll2007l) magnetic O 



star HD 191612, using an inclination an gle i = 50° and = 50°, consistent with the 
magnetic geometry constraints derived by Wade et all (T201ll) . Fig.[5]shows that both the 



peak flux at phase and the extended minimum around phase 0.5 are well reproduced. 
The time modulation is caused by the observer's changing projected surface area of 
overdense Ha emitting wind plasma as the star rotates. At maximum emission (phase 
0), the DM is viewed essentially from above the magnetic pole, whereas at minimum 
(phase 0.5) it is viewed almost edge-on. Thus the DM in the latter phase has a smaller 
projected surface area, leading to weaker emission. 
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The good overall agreement between observations and simulations obtained for 
HD 191612 is encouraging, and provides strong support for the basic concept of the 
dynamical magnetosphere model. We are currently in the process of extending the 
modeling presented here to more magnetic O stars with /? + < Ra < Rk that also show 
periodic Balmer emission, like HD 57682 (Grunhut et al., submitted to MNRAS) and 
6 l Ori C. 

Acknowledgments. J.O.S. gratefully acknowledges funding from NASA ATP grant 
NNX11AC40G 
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